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Abstract 


Introduction 


Steady-state solutions of the Reynolds-averaged 
Navier-Stokes (RANS) equations were computed using 
the Cobalt flow solver for a constant-section, 
rectangular wing based on an extruded two-dimensional 
glaze ice shape. The one equation Spalart-AUmaras 
turbulence model was used. The results were compared 
with data obtained from a recent wind tunnel test. 
Computed results indicate that the steady RANS 
solutions do not accurately capture the recirculating 
region downstream of the ice accretion, even after a 
mesh refinement. The resulting predicted reattachment 
is farther downstream than indicated by the 
experimental data. Additionally, the solutions computed 
on a relatively coarse baseline mesh had detailed flow 
characteristics that were different from those computed 
on a refined mesh. However, the results predicted using 
the baseline mesh were conservative in the sense that 
they indicated a more severe reduction in lift than the 
results computed on the refined mesh or the 
experimental data. Steady RANS solutions were also 
computed to investigate the effects of spanwise 
variation in the ice shape. The spanwise variation was 
obtained via a blending function that merged the ice 
shape with the clean wing using a sinusoidal spanwise 
variation. For these configurations, the results predicted 
for the extruded shape provided conservative estimates 
for the performance degradation of the wing. 
Additionally, the spanwise variation in the ice shape 
and the resulting differences in the flow field details did 
not significantly change the location of the primary 
reattachment. 
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There are two distinct applications related to the aircraft 
icing problem in which computational fluid dynamics 
(CFD) can play a significant role. 1 First is the 
prediction of ice accretion using software such as 
LEWICE 2 coupled with a viscous flow solver. The 
second application is a detailed flow field analysis to 
determine the effects of ice accretion on aircraft 
performance. We focus on the second application, 
“icing effects,” in this paper. 

Recently, numerous simulations of the flow fields 
associated with airfoils and wings with ice accretions 
have been reported for “icing effects” studies. Chung, 
et. al, 3 performed a CFD analysis of the ice 
contaminated wing surfaces of a turbo-prop aircraft. 
Their analysis was performed to obtain qualitative 
trends and to provide insight into the aerodynamics that 
may have led to a loss of control of the aircraft. Dunn 
and Loth, 4 Dunn, et. al, 5 Kumar and Loth, 6 and Pan, et. 
al, 7 computed the effects of simulated spanwise ice 
shapes on airfoil aerodynamics. Their simulations 
showed that for two-dimensional computations, detailed 
flow characteristics such as pressure distributions and 
reattachment lengths showed reasonable correlation 
with experimental data, as did integrated quantities such 
as lift, drag, and moment. Notable exceptions to this 
trend were that the agreement was degraded once the 
flow on the upper surface was fully separated and that 
the maximum lift values were not well predicted. In all 
cases mentioned above, the steady Reynolds averaged 
Navier-Stokes (RANS) equations were solved to 
compute the flow field. Kumar and Loth 8 performed 
two- and three-dimensional detached eddy simulations 
(DES) for iced airfoil/wing configurations. The 
agreement of time-averaged results with experimental 
data showed some improvement when compared to 
steady RANS computations. 

Although computational fluid dynamics tools have 
made great strides in recent years, there are several 
unanswered questions regarding their capabilities for 
predicting the flow around wings with ice accretions. 
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The ability of geometry/mesh generation tools to 
adequately address the unique needs of an iced wing 
flow simulation is also much in question. There are 
significant issues related to the iced wing surface 
modeling including questions regarding the quality of a 
nonuniform rational b-spline (NURBS) representation 
of the surface. Additionally, within the constraints of 
accuracy, efficiency, and flexibility, questions remain 
as to what type of mesh is appropriate for a complex 
iced-wing configuration. Finally, there are numerous 
questions related to turbulence modeling, e.g., regions 
of massive separations and surface roughness. 

The effort reported here had two objectives: 1) to 
evaluate the effectiveness of current computational 
fluid dynamics technology for predicting flow fields for 
a wing with a glaze ice accretion and 2) to investigate 
the effects of spanwise variation of ice shapes on the 
resulting computed flow fields. Steady-state solutions 
of the Reynolds-averaged Navier-Stokes (RANS) 
equations were computed using the Cobalt flow solver 9 
for a constant-section, rectangular wing with an 
extruded two-dimensional glaze ice shape. The one- 
equation Spalart-Allmaras turbulence model 10 was 
used. The results were compared with data obtained 
from a recent wind tunnel test. 11-12 Two-dimensional 
predictions made for the same airfoil/ice accretion are 
presented in a second paper in this session. 13 

The basic ice shape considered was the 22.5-minute 
glaze ice accretion on a GLC305 airfoil, which is 
denoted as the two-dimensional 944 glaze ice shape. 14 
The effects of spanwise variation in the ice shape on the 
resulting flow field were also investigated by 
computing flow fields for synthetic ice shapes that were 
generated by merging the two-dimensional ice shape 
with the clean wing using a sinusoidal blending 
function and comparing these results with clean wing 
and extruded wing predictions. 


• Specify the artificial boundaries needed to 
define the computational domain, e.g., the 
outer boundary and the side boundaries 

• Generate a mesh and specify boundary 
conditions on the bounding surfaces of the 
computational domain 

• Generate a mesh in the interior of the 
computational domain 

• Generate a flow solution 


Since the wing employed for the test program was 
mounted between walls, the geometry modeled in this 
study did not include wing tips. 

Geometry Modeling 

We developed a technique to generate a wing with a 
synthetic three-dimensional ice accretion given an 
airfoil definition and a two-dimensional ice shape. The 
algorithm interpolates between the clean shape and the 
ice shape using a sinusoidal weighting that depends on 
the spanwise location. 


We now describe the basic algorithm. Given a point on 
an airfoil with an ice accretion (xi cei i,y ice d) and the point 
on the corresponding clean airfoil (. x c i ean ,y c i ea n ) 
definition that is closest to {x i(;ed ,y lced ), interpolate 
between the iced shape and the clean shape using a 
factor that depends on the spanwise position 


= x clea n + i X iced ~ X rlean) X ^(z) 
y new = y clean + ()W “ >«*«») X <t(z) 


(1) 


where o(z) is defined as 


cr(z) = max] 0.0,1 + — 5 cos[2 (on 


2 V V 



(2) 


We first provide a description of the computational 
problem and present the method used to generate the 
synthetic ice shapes. We then briefly describe the mesh 
generation process and the flow solver. Finally, we 
present the results obtained to date. 

Problem Statement 

The prediction of the aerodynamic characteristics of an 
iced wing is a complex problem that involves several 
steps. In order to motivate our approach, we now 
provide a problem definition. Given a geometric 
definition of a wing with an ice accretion: 

• Develop a representation of the surface that is 
suitable for generating a surface mesh 


8 is the deviation from the clean surface, z is the 
spanwise position, co is the frequency of the oscillation, 
and b is the total wingspan. Note that the max function 
ensures that o(z ) is nonnegative. The function 
multiplying the deviation <5 varies between 0 and -2 so 
that the maximum value of o(z ) is 1 and the minimum 
value is 1 -8 with a sinusoidal variation in between. If 
o(z) is unity, the iced shape is returned. If o(z ) is zero, 
the clean shape is returned. It should be noted that, for 
complex ice shapes, a ray extended from the clean 
airfoil surface in the direction of the local surface 
normal may intersect the ice shape more than once. The 
above procedure may produce a "folded" surface in this 
case if larger values of the deviation parameter 8 are 
employed. The 944 glaze ice shape is one such shape. 
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Figure 1 shows four geometries. The “4-cycle” and “8- 
cycle” cases were generated using 3=0.2. These cases 
were selected because the resulting ice shapes exhibit 
spanwise variations that resemble realistic three- 
dimensional ice accretions. Figure 2 shows the effects 
of varying 3. Figures 2(a) and 2(b) show valid surfaces 
generated using 3= 0.2 and <5=0.5, respectively. Figures 
2(c) and 2(d) show invalid surfaces that are generated 
using (5=1.0 and 3= 1.25, respectively, and are included 
only to show the effects of 3 variation. 

Mesh Generation 

For each geometry, a structured surface mesh was 
extruded from the section description according to Eq. 
(1). This surface mesh was then converted to a NURBS 
representation using the mesh generation software 
GUM-B. 15 The baseline (coarse) surface and volume 
meshes were generated using SolidMesh. SolidMesh is 
an interface to the unstructured surface and volume 
mesh generation software AFLR2 and AFLR3. 16 
AFLR3 uses an advancing front algorithm to insert a 
point in the mesh. The point insertion is followed by a 
local reconnection to improve mesh quality. The refined 
surface and volume meshes were generated using 
GridTool 17 and VGrid. 18 VGrid uses an advancing front 
algorithm to generate a tetrahedral volume mesh. For 
the all-tetrahedral meshes generated using VGrid, the 
near-body elements are converted into prisms 
producing a mixed element hybrid mesh. SolidMesh 
produces a hybrid mesh with prisms and tetrahedra 
automatically. These hybrid meshes were employed 
because of their potential for improved efficiency and 
accuracy in comparison to unstructured tetrahedral 
meshes. The GridTool/VGrid combination gave 
somewhat more user control of the point spacing than 
SolidMesh for the wing surface as well as in the flow 
field. 

Table 1 shows statistics for the five meshes employed 
to generate the results reported here. In all cases, the 
distance to the first point off the wall was defined so 
that an average y + value less than 0.1 was obtained. 
This value is well within the recommended values for 
the turbulence model employed. 

Figure 3 shows cross-sections through the baseline and 
refined meshes generated for the wing with the 
extruded ice shape. Note that these are cutting planes so 
that the line segments that appear represent 
intersections of faces with the cutting plane. However, 
the connectivity of the mesh is evident and the prism 
layer near the body surface is clearly visible. Also, 
Figure 3(a) shows the faceted surface represented by 
the triangular surface mesh on the extruded wing. 
Notice that even though the surface definition is two- 


dimensional, the resulting faceted surface exhibits a 
spanwise variation due to the unstructured nature of the 
mesh. Additionally, due to the manner in which the 
spacing is controlled in SolidMesh, there is a spanwise 
variation in the spacing. The volume mesh in the region 
downstream of the upper horn is relatively coarse. 
Refining the surface mesh reduces the faceting and 
improves the mesh density in the region where flow 
separation is anticipated as shown in Figure 3(b). 
Additionally, a more uniform distribution of edge 
lengths is obtained over the surface using the 
GridTool/VGrid combination. 



# Nodes 

# Faces 

# Cells 

Clean wing 

630,285 

4,711,884 

2,150,999 

Extruded wing 
(baseline) 

890,958 

5,299,738 

2,258,447 

Extruded wing 
(refined) 

1,828,711 

17,327,367 

8,341,019 

4-Cycle wing 
(£=0.2) 

904,680 

5,368,010 

2,285,734 

8-Cycle wing 
(£=0.2) 

903,654 

5,360,055 

2,281,964 


Table 1. Statistics for meshes employed in this effort. 


Flow Solution 

The flow solver employed in this effort is the 
noncommercial version of Cobalt. 9 The Cobalt flow 
solver was designed for general unstructured meshes. It 
employs a nonlinear Riemann solver for the inviscid 
flux computations and can be run either in explicit or 
implicit mode. In all cases reported here, the solution to 
the RANS equations was obtained using first-order, 
implicit local time stepping and does not in any way 
represent a time-accurate solution. Second-order spatial 
accuracy was obtained using a linear least-squares 
reconstruction of the data. Several turbulence models 
are available including the Spalart-Almaras one- 
equation model 10 that was used for the computations 
reported here. A slip boundary condition was applied 
on the artificial side boundaries. 

Results 

In this section we present a comparison of predicted 
steady RANS results with experimental data 11 ' 12 for the 
rectangular business jet wing with the extruded 944 
glaze ice shape. We also include a brief discussion on 
the effects of spanwise variation of the ice shape. 

We computed steady RANS solutions using the Cobalt 
flow solver with the Spalart-Allmaras turbulence model 
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Figure I fa). CEean GLC305 wing settioia 



FEgure EslnJiied U44 j)*ze ice wing section 



Figure ] {d)_ " S-cydc 7 ’ wing section 
S=Q2, £j=S 

Figure I . Iced wing «H*flgura1ions <kmoustrt*fog (he 
effiirts of tv variation 



Figure 1 (a). ’'S-cycle'' wing section 
’ S=OX 07 =S 



Figure 2(h). "S-cyde T ‘ wing suction 
S=0 .5, ta = 8 



Figure 2(c). "i-Gycie"' wing section 
£*=L A *»=« 



Figure 2(d)- “S-cycle" wing: stslbn 

^1.25. w=8 


Figure 2 . Effect of d’ variation on resuliing wing shape, 
Figures 2(c) end 2(d) represent invalid geometries due 
to surface “folding. 
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Figure j{a). Basslira mesh for extruded Lee shape 


Figure ZhTd). Refined mesh few extruded ice shape 


Figure 3. Mesh in cutting plana fiw extruded 944 glaze 
ioe sFiaps - i» baseline mesh and (b) retimed mesh 

for die following, conditions-- M=0,l2 h RealXl.Sxl 0 6, 'in 
which., with a chord Length of 0.9144m, yields 
Re=3.5j(]0 6 1 and angles of aua-;);. of 0f\ 2", 4* h and d D . 
The clean GLC3Q5 wire;; cases correspond to- Run 9 in 
the experimental dato, ‘Hie extruded ?44 ice shape casss 
correspond to Run 4 1 . 


shows the the L? norm of the residual. The nosiduaL 
decreases slightly more than oire under of fita^nirude 
during the computation. Figures 4(h) amd 4(c) show the 
lift *dd uxjal Ibree histones, respectively. Notice fftar in 
both casts, the solution appears to have reached a 
steady state. 


Figure 5 shows compurisorjs ofccrnpuled result* w-ifo 
experimental data. Figure 3(a) shows a comparison of 
the predicted and experimental Liti coefficients obtaiueri 
hy integrating the surface pressures. In this ease, there 
as excellent agreement between ihe predicted results 
and the ex pe rim ear a I dlara. This is to be expected for 
these rcln-tivcty law angle of attack cases, Figure 5(b) 
shows a comparison of predicted and ox peri mental drag 
curves. The drag coefficient was obtained for the 
prediction by integrating the pressure and foe viscous 
stresses over the surface of the airfoil The experimental 
drag force coefficient waj obtained by integrating a 
w.sl:c nuvey. In this ease, the data docs not show 
agreement similar to *haa obtained for the lift results. 
This occurs because. For (he streamlined body under 
cMi*id4re4cm # lew wigle of attack, the d«g is 
primarily due to viscous. stresses. Finally. Figure 5(c) 
shows a comparison of lbs predicted and expert mental 
pressure coefficients at am angle of atlacit of 6". As 
suggested by the lift curves, the agreenvjnt between the 
prediction amd experiment is very good. There wjeie Little 
spanwi&c variaiion in the computed jorfect pressure, 

Fxhudiad Wing: Coifina-H.^n with Exp^rimriilal fJala 


Ail cases reported below were run on £4 proeossore on 
the EMPIRE cluster at the ERC at Mississippi Slats 
University. The EMPIRE cluster is a supercomputer- 
class clusteT of workstations consisting of 105® one- 
GHz or better Pentium III preeesflots eseh with one or 
more 06 Of RAM. Solutions computed -on this baseline 
mesh for foe extruded wing required approximately i 
CPU hours per thousand time steps. Solutions 
computed on the refined mesh foe foe -jxtmdisd wing 
OOtffigMraiion t«i aired approximately 1 b CPU hours per 
thousand time steps. 

Clean. Wing: Comparison with ExuLiimoma! Data 

For the purposes of establishing a baseline level of 
confidence for the methodology employed here, we first 
compare predicted result* with exp*rim»rtta.l dirt for 
the cLeUn wing. i.e.,. a rectangular wing with the 
GLC305 section and sn aspect ratio of unity, 

The first point to consider is- the- oonvwgofice of the 
soliffion. Figure 4 shows several convergence measures 
for the 6° angle of atrnck imputation- 6 is 
representative of the other compulations. Figure 4(a) 


We now- compare result* ctfmpitted for foe extruded 
wing based on the two-dimensional $44 glaze i-ce 
shape. 1 " Ccmtpuwions. were perforated on ihe meshes 
denoted as ^'extruded wing (baseline)" and "extruded 
wing (retimed)" in Table 1 . 

Figure 6 shows. (he convergence histories for she 
extruded wing computation for an angle of attack of 6° 
using foe baseline flash. The Lj residual, foe lift forte, 
and the axial force reached asymptotic values as shown 
in Figures 6(a), 6(b), *ftd b(c), respectively. This 
computation is representative of the computations for 
other angles of attack using this me&ri 

Figure 7 shows the convergence histories the 
extruded wing computation for an angle of attack, of E 5, 
using Ac refined mesh- Figure 7m) shows that the L? 
residual has reached an asypmtotic value. The lift force 
and axial force Kivu also reaped .asymptotic values as 
shown in Figures 7(b) and 7(e), respectively,, This 
convergence history is representative of the 
computations for other angles of attack using this mesh. 
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Figirre 5(a). Clean wing - lift coefficient vs. angle of 
attack 


4ta- 



Figure 4 (b). Clean, wing li ft history 



F^pw* 3(b)- Clean Wing - dr.ig coefficient vs. angJeof 
at tael; 



FtfiOn* 4(c), Clean wing - aflial foie e history 



Figim: 5(c). C l'.'u.-i * jfig ctumiwise di stri butiwi of 
prssure coefficient - fltr6* (midspan) 


Figure 4 . ConvE^gence of clean CSLC3fl5 'wiflg 
conipiiaticin - oefS 11 


Figure 5 Comparison of predicted risi.il Is with 
experimental daia few the e xirwfed GLC 305 wing 
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! IU.il ;■ t>! 

hurt iiiwi 

figUrt Estrtded wing - Li residua] hisiocy 
(baseline mesh.) 



(refined mesh) 



(baseline mesh) 



Figure Extruded wing - wing axial force history 
(baseline mesh) 



(refined mesh) 



(refined mesh) 


Figure -6. Convert edl££ of extruded $44 glfcz e ice shape Figure 7. Convergence of extruded 944 glnze Lee shape 

eortHJuMfcion - Ot=h t (baseEbe mesh) compulation - if=C b (refined mesh) 
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Figure 6 shows a comparison of Dw wing lift md drag 
coefficients pradieacd using the base tine mesh and She 
re lined mesh with experimental data ft* ait angle of 
attack range of CT t* fi° in 2" increments. As shown in 
Figure S(a), Ihe lift coefficient ai {F is well predicted. 
Howevex, a* the angle of attack is increased, the I ill is 
Wlderp re dieted relative la the experimental data. 
Addiitiofially, die break in Ihe predicted lift curve 
s'ope is eharaclerbtic of the ‘ near stall 1 * behavior of the 
wing suggesting thin the predictions are overestimating 
ihe extent of (he flow separation. Oh tire raher bond. as 
shown in Figure 8fb) h ihe predicted drag cocffkiem 
shows very good agreement with Die experimental dala. 
This agreement be attributed to the fact that the 
drag ia primarily oomprx&ed of rom or pressure drag for 
this ccmJljUrtfion i e., the ice shape is a bluff body. 
Howev er, the vary §<>04 agreement between the results 
predicted using the baseline mesh and die es; peri two I4l 
dala appears ta he fortuj-uxis The agreement is not as 
,3}f>od f®T predictions made using the refined mesh. h 
should be noted lTielv for each else considered, there are 
only very 1 slight spaowise variations in the pressure 
distribution. 

Figure 9 shows a comparison of the predicted midspan 
pressure weflfltiente id: the baseline mesh and ihe 
refined mesh end experiitwOLiD dala for an angle of 
attack range of 0“ to £> n in 2* increment. Few the most 
part, the tower surface pressures are well predicted 
nsing eitiier the baseline mesh or the refined mesh. The 
differences in the predicted upper surface pressure 
eoeffKiwus are consistent with those in ihe liff 
coefficients. Mesh refinemem improves agneemeot with 
ihe experimental data somewhat and Indicates a 
strengthened recirculation region. These results 
Suggests Drat further mesh raftneCKGl may be 
necessary. However, it may be lhal a steady RANK 
solution w ill not accurately tapujre tht T 2 ciTi;uUiing 
region and an alternative procedure such as an unsteady 
ffeLaiched Eddy simutHiion (DES) 1 * may be necessary 
for this case. These remits are consistent with those 
reported by Chi, et, al. IJ 

Vl r e now include a detailed comparison of computed 
flow field velocities, (the u- component of the velocity) 
wall experimental dHta ohrained using; [he sp-liL-film 
technique as reported by Bruerea, et al. Hie iqriit-film 
technique is a variation of hot-wire ahemometry. 
Metallic films are deposited on apposite sides of a 
small cylinder. The cylinder is then placed 
pemcr.di euiar to the primary flow direction of interest 
and heated using an electric current. By monitoring the 
vo]r ages needed to maintain constant current, the 
magmMde and direction of flow jp the primary 
direction, u, can be determined. The magnitude of the 
v-component of velocity may also be ascertained, but 


not its direction. We do nat include any coin pa crises 
with flit v -component and. evocc our results are steady 
RANS oomputations, we consider .comparisons only 
time -avenged mean flow values. 

Figures 10-12 show comparisons of the predict <sd u- 
componefit of velocity v. iih experimental dam fi>i- bdll 
she baseline mesh snd the refined mesh for angles of 
attack of tf 1 , 4°, and 6 11 . ]t should be noted that the- Color 
scales employed to display the predicted data and the 
experimental dala are somewhat different- Using these 
images, the reattacbmcnC location of the primary upper 
surface flow separation may be cttcratHed by legating 
the po&itxm on the chord at which the contour u=U 
rotcrsects the upper surface uf Lbe wmg, Th? general 
trend exhibited m these images ts tliai, at an angle of 
attack of (F h ihe rcMfedtident position is well presided 
using cither the baseline mesh or the refilled Mnesh. 
However, as Die angle of attack increases, the raJUrait of 
the separated region is overpre dieted until the fl<sw is 
completely separated on the upper surface at an angle 
of attack of 6". Mesh refbienrwiii does impeov-c Ihe 
prediction relative to the baseline mesh. How ever t,, the 
extern of the separated region is still overesti mated!!. The 
™a dissipation present in the baseline mesh solul ion is 
evident from the increased spreading; of Die shuar layer 
in comparison to results predicted using the refined 
mesh. However, Dais spreading is not at ttven as 
expected given the relative coarseness of the baseline 
mesh. The seemingly anomalous velocity wntours thar 
appear in the experimental dala Just downstream *txf the 
horn are artifacts from Die pr ww employed to 
generato the cotitora plots. 

Figures Id und ]d show predictions of separation, (red) 
and attachment (bluo) lotatijOTfe for #0 angle of ^tuid; 
range of if ic 6 C in 2' 1 increments for the baseline mesh 
and the renoed mesh, resp^tively. The te chut iq lie of 
K«lwrighL i ‘'' is employed to locate surface mesh 
eleraents m whieh a ponjaiial -nifraniijoii or ilHtfrinwit 
Line -ctas^tt. No aUempI is miide here to reconsirtf ct Ihe 
Hctuai cur r es. KenwTight’s todinique ii ba^d oei a 
phase plane analysis of critrc-al points in the loeity 
field ^ciose 15 ' to the Arfaue aiid h a® suth. is sabi to 
aaomalous results ifthere is “Ttoise” in tbe data. 

The refiired mesh results show fairly w-cll rcsoh'erdl flow 
features iocludipg &$co*idary chordwlAc KpftrftiCTU and 
*OfiWpohdmg atHchmenls clearly indicated just down 
stream of the bom (ihe closely 5isaced h roughly parallel 
blue apd red curves Located jus aft of the leading esdge). 
The refined mesh feso-o. show a relatively tou ifpnii 
primary reaLLJchnncnl across the span in each case. The 
baseline mesh dots hot rssclvc this secn> ttdary 
separoinn at any angle of attack. Additionally, the 
iocalion of tha pHtHSry reoitaichraenl predicted Using 
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Figure S(a). Extruded wfeg lift coefficient vs. angle of Pignut S|b>_ Extruded wing drag coeffldettf vs, angle 
afflaclt of attack 

Figure 3. Comparison of predicted fierce Coefficients whb experimental data fertile extruded! ifcid gla^e its shape 

(baseline mesh and rersned mesh) 




■ — ■ I L'. 7 « ■ —r/> J.:-, 

«~l ! : FT-*- 1’ ■ I W«i ■■ V I 

HF^V iinfagJ irb 1 ' 


*k 

Figure (*(h). Extruded wing pressure distri sution — ct=2 


-L— ■ 1 1 I I 

W OlA □.*■ QJ 

ifc 


Figure ■'■>!». Extruded wing pressure distribution - a=ff 


0 



2 1 ! ■ I 1 1 i 

a HI 0>I i* S>JS 1 

ah 

Figure 9(c). Extruded wing pressure distribution - c=^ 


Figure 9 (d), Extruded wing pressure distribution - 01=0" 


Figure 9 , Comparison of predicted pressure oocfficiesiis {midspan) with experimentsE data for Hite extruded 944 

jliM ice shnpe (baselme mesh and refined mesh} 
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xJSn 

Figyjc 10 (c). Extruded ’Air.g (cKpcjimcutal data) -e-velocity wntotn^ 

Figure 10, Comparison of-' predicted e-velodEy concurs (tnidspan) with experimental <:jta ri>r l^h excrud&d 944 glue 

ice shape fbascLine mesh and refined me ah) - sj =alV J 
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Figure 11(b). Extruded, wing (baseline mesh) - u-vdodly cmHoujs 



Figure L l|Jb}. Extruded wing (tefli-ed mad i) - u-velocity cuntowt 



jrft 

Figure 1 1(c). Extruded wing (Btperimeftttl date) -u-vrlodty contours 


flours 11. Comparison -of predicted u-ve bd(y 4 ertiroure (tnddspaa) ^ilh experimental data for the extruded 
glaze ic* shape (haw | jm nncsb and refined mesh) - Qt=4 p 
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Figure L2[a). Extruded wing (ba«Hne mesh) - u-veEocity contours 



Figure ]2(b). Extruded wing (refitted mesh) - u-Ye Booty contours 



a>c 

Figure 12 (e). Extruded wing (experimental data) — u ^velocity contours 


Figure 12. Comparison of predicted u- velocity contours i mid-span i with experi menial dalu for the extruded 
9A4 glaze ice strips Cbasoline mesl and refined mesh) - 0 = 6 ’' 
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Figure S3fc). Upper surface separalion and rwBachjnaUt lines (left) and streamline treecs | right) 



Figure 13(d)- Upper surface separation andreattaclunent lines (left) and streamline Iraces (fight) - <x=tf* 

Figure 13. Predicted upper surface separation {red) and attachment (Hue) linen and sLream'ii e lt^ccs colored fry 
the u-component of velocity - extruded 9*4 glare ice shape (frasclir.e mesh) 
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Figure 14(b). Upper surfat*- wparelkrti -and reauadimKH liacs (left) «nd streamline e races (right) - ct~2 u 




Figure 14. Predicted surifce* seporrunon (red) and attacbmait (blue) ling and streamline traces rolored by 
tbe u-ccTT.poiiir-t of velocity - lotfvud&d 944 glaze ice shape (refined mtsb) 
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doss not exhibit the same ufofOTtmy that is present in 
the refined ffltcsh computations me^t for l he 6 y case 
which is fully separated. Wt surges; (hen (he coarser, 
nonimifOTtfl surface tnesh employed in the ha^iinc 
mesh {see Figure 3fs)} is reasonable for the apparent 
ttomLTJ tOTmi ly of the baseline mesh flow fields. 
CMicLuskma drawn rajja«ting (be attachment locations 
are consistent with those drawn feom the velocity 
contours for foe 0^, A\ and etscs. 

Also shown in Figures 1J and 14 are instantaneous 
streamlines tit eiich: angle of attack for both meshes, The 
nonuniformity in foe baseline mesh flow field is also 
evident ip the stream lines. At each angle of atLacV, the 
streamline traces for the baseline mesh solution show 
considerably more three- dim.cn signal behavior tbiin Ihe 
streamline traces for the refined mesh solution. 
However, we do not suggest that the baseline mesh 
more accurately models an teed wing flow field. The 
random faceting, which occurs because of the 
orientation of die surface triangles, appears at the 
highest frequency represefiiablo by the mesh. Allhough 
the faceting of the surface shown in Figure 3(a) docs 
more dosely resemble an iced wing, surface than (be 
smoother description, shown in Figure 5(b), the flow 
features associated with this faceting would atrt be 
resolved on this, mesh K is difficult to state 
conclusively that the chord wise curves represent real 
fcamres in ihe data. It is likely that they represent 
inadequately resolved secondary [low- features. 
Additionally, the predicted results show that once the 
flow becomes folly separated, these seeondaty 
structures lose tfiii r coherence. 

It should be rimed that, although the details of ihe ifow 
field were not well predicted using Hie baseline mesh, 
the predictions ;q?peuT to be conservative. That ia, the 
lift degradation predicted by rb; baseline mesh 
solution* is more severe than the experimental da*a. 
The is primarily due to the more extensive flow 
separation predicted using (he baseline mesh. 

Effects of Swtwise Varintiou 

We now consider Hie effects of spanwisc variation in 
the -ce shape. The configurations employed here are Ihe 
four eon figurations showxi in Figure I - foe clean 
GLC305 wing, foe extruded $44 glaze ice shape, the 
"i-cycle" wing,, and (he "8-cycle 1 ’ wing. The 4-cycle 
and 8-cyde cases with S=CL2 represent reasonable 
three-dimensional ice shapes. The meshes for these 
rases* described in Table 1, arc roughly equivalent, in 
resolution, to the ha** K m mesh for (he extruded wing. 



Figure 15(a). Wing lift eocJYteienl 



U.l jT .•: J LA. l I Jfj,l 


Figure 15(b). Wing drag coefficient 

Figujc 15. Comparison of predicted wing lift and drag 
coefficients for four different configurations clean 
GLC3Q5 wing, extruded $44 glaze Ice shape (baseline 
mesh}, 4-eyde wing with nnd 8 -cycle wing with 

$=cn 

figure 15 shows o comparison of the lift and drag 
coefficients predicted lor each of the four 

Configurations across an angle of attack iaiigi of 0 D to 
in 2° increments. Of x\ j,ni ficunce here is the iaci that 
the clean wing and extruded wing hound holh the lift 
aud drag behavior of foe other configurations. This 
implies lhat the extruded wing provides a cgnjervative 
estimate for the predicted l ift and drag. On the surface, 
(his result appears reasonable since the spanwisc 
variation in the ice shape effectively lessens foe volume 
of foe accreted ice and produces a smaller blockage 
effect. However 1 , this conclusion has not been 

substantiated for other configurations and should not htt 
extrapolated 3i this iime. 
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FlgLiTCS 1 6 EJld 1 7 sl»w predictions of SCpaiBtLOD (red) 
and attachment (bin:) locations for an angle of attack 
range of 0° co in iticreftwdis for tbs 4-cycle wing 
and the fi-cycle wing, respectively. in general, the 
results indicate the presence of secondary stimetuncs 
lhat are Introduced bv the s^anwi^- «ria:?cni in the Lae 
shape. The structures, interact wish the Elow field to 
produce an irregular pri ns ty reattach isnem region. Not 
surprisingly,. til ere is i tfuaai -periodic variation that 
appears to be associated with the span 1 , civ? viri at ton of 
the ice shape. [1 should be noted, that these secondary 
structures are not adequately iesob^4oD these melws. 
Also included ia Figures It and E7 arc the 
instantaneous streamlines foe each u-i£]e of attack. Ae 
cxptL-Ui-J, the streamlines show considerable three- 
dimensional flow. The streamlines mdicule the presence 
of wra p I ex flow pu items in the region just duvvais.rca-in 
of the hom, even at thfc lower ingles, of attack. 
However, ike spao^tse flow’ docs not appear io 
significantly alter the global characlenstics of the flow 
field to particular, althougli the overall shop* of die 
primacy reattach men! varies between Figures 13 
(extruded 944 glaze ice shape),. 16 (4 -c vela w ing), and 
17 f S-cyde wng), its locution is, essentially unaltered. 

Conclusions 

The results computed as pad of this effort indicate that 
tlw steady RAKS solutions do not accurately capture 
the recirculating region downstream of Ihe 944 glaze 
ice accretion, eve® after a mesh refinement, and that the 
resulting predicat'd reattecbfficiil is fa diner downstream 
thao thtr; indented by the experimental data. 
Additionally,, die detailed flow characteristics m the 
•solutions computed on the coarse, baseline itie^h were 
different than those compuled on tbs refined, mesh. It 
should be noted lost the results :>iestri:us£ here do not 
really represent a grid sensitivity study since grid 
independence of the solution has itot been 
demonstrated. However, white Cbe detailed flow 
characieri slics on the two meshes were different* the 
wing lift values far flic two meshes are quite similar. 
Further, the results predicted using the btethtit imh 
were conservative- in the sense that they indicated a 
more severe reduction in Sift than the espcTirnenml d : ita, 

A spnnwise variation to ice shape was obuined by 
developing a hletldHng function that merged the 944 
glttzc ice shape with the GLC3QS clean section using ii 
sinusoidal span wise variation. Limited results indicate 
■h:it. :br this case, the results predicted by the extruded 
shape provide cohscrvati ue estimates lor the 
performance degradation of the wing. Additionally* die 
spanwise variation in the ice shape and the resulting 
differences in the flow field devils did not significantly 
change the: location of the primary' re-attachment. Tfiesse 


results suggest Hie influence of moderate spnnwise 
variation to (he tcc accretion may be only minor for 
gross quantities sucti as Lift and drag 

The results reported in (his paper indicate that (he 
agreement between cmviputed result* and eXperiivneflial 
data is degraded for cases with large regions of upper 
surface flow separation and are consistent with those 
reported for two-dimensional computations:. 1 " "These 
Jesuits suggest til at funher refinement may he needed li- 
Ihal an alternative approach, such as unsteady DE5, 
may be necessary and liiai further study of thasc issues, 
is warranted, 
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Figure ] 7(b), Upper sunset separaiian and rtattaebrurnt lines (left) *ud streamline traces (right) - u.=l z ' 
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figure i 7(d). Upper surface separaibn and niaitachn^nt Hrte=& (left) and Efteajnlrae traces (right) - a=ff* 

Figure 1 7. Fredicted upper surfajes separation (Vwl)und anachment (blue) ibes and streamline Paces colored by 

the U'tUfaporiEii ofirekKity - S-c>cLe wing with 3=0.2 
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